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Exact and asymptotic formulas relating to dynamical correlations for over- 
damped Brownian motion are obtained. These formulas include a generaliza- 
tion of the /-sum rule from the theory of quantum fluids, a formula relating the 
static current- current correlation to the static density-density correlation, and 
an asymptotic formula for the small- behaviour of the dynamical structure 
factor. Known exact evaluations of the dynamical density-density correlation 
in some special models are used to illustrate and test the formulas. 



1 Introduction 

A number of recent papers [1-4] have considered various aspects of the density and cur- 
rent correlations in the Dyson Brownian motion model [5]. This model refers to the 
overdamped Brownian dynamics of the one-dimensional one-component log-gas. When 
the initial state is given by the equilibrium state, the model is equivalent [1] to the 
ground state dynamics of the quantum many body system with 1/r^ two-body interac- 
tions (Calogero-Sutherland model). One motivating factor for the current interest is that 
the Dyson Brownian motion model specifies the eigenvalue probability density function 
(p.d.f.) for certain ensembles of parameter dependent random matrices (see e.g. [6]). 
These random matrix ensembles are used to study the parametric dependence (typically 
as a function of magnetic field strength) of quantum energy spectra and fiuctation phe- 
nomena in quantum transport problems [7,8]. Another motivating factor is that there 
are now a number of exact results available for the dynamical density-density correlation 
[4,7-9]. 

In this paper we seek to extend the considerations of earlier works by considering 
properties of dynamical correlations for overdamped Brownian motion described by the 
general Fokker-Planck equation 

dp ^ , ^ ^ d fdW ^ , d \ 

lTf = Cp where £ = ^ — -- + /^-i— . 1.1 
or jr^ oxj \ oxj oxj J 
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The Dyson Brownian motion model is specified by this equation with the particular 
potential 

1 ^ 

W = -J2^^j- E ^og\xk-Xj\, (1.2a) 

j=l l<j<k<N 

or its periodic version 

\<3<k<N 

We begin Section 2 by reviewing properties of (1.1). We then clarify the meaning of 
the current (notice that (1.1) is independent of the velocities). Next a generalization of 
the /-sum rule for quantum fluids [10] is given, as is a generalization of the formula of 
Taniguchi et al. [3] relating the initial density-density and current-current correlations. 
In the final two subsections of Section 2 we consider the hydrodynamic approximation to 
the microscopic density fluctuations, and derive from it the small-A; asymptotics of the 
density-density and current-current correlations, which in turn imply certain sum rules 
and further asymptotic formulas. The analytic formulas of Section 2 are illustrated and 
tested on some exact results for the density-density correlation in the Dyson Brownian 
motion model (i.e. (1.1) with potential (1.2)), and for (1.1) with potentials closely related 
to (1.2). Our results are summarized in Section 4. 



2 The density and current correlations 

The Fokker-Planck equation (1.1) relates to the overdamped Brownian motion dynamics of 
a classical gas with potential energy W in contact with a heat bath at inverse temperature 
/3. The function p is the p.d.f. for the event that the particles are at positions Xi, . . . , xjv 
after time r. As written in (1.1), the particles are assumed to be conflned to a line, 
however the same equation applies in higher space dimensions if we simply replace Xj by 
x^f^ (the components of the position coordinate Xj) and assume summation over a. 

The dynamics specified by (1.1) is known [11] to be equivalent to the dynamics specified 
by the coupled Langevin equations 

^^ = -15-+^^M = l....,iV) (2.1a) 

where the random force J-'j (r) is a Gaussian random variable with zero mean and variance 
given by 

W:F,{t') = ^S,,S{t-t') (2.16) 

(the average denoted by an overline is with respect to the random force). In particular 
this means that the correlation functions can be calculated using either (1.1) or (2.1). 

In the Fokker-Planck formalism the correlation functions can be specified in terms of 
the Green function G{xi'\ . . . ,x^^^\xi, . . . ,xn]t), which is by definition the solution of 
(1.1) subject to the initial condition 

p{xi, ...,XN;r) = n S{xj - xf). (2.2) 
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Thus for observables A-,- and Br (e.g. At — nr{x) := J2jLi3{x ~ ^ji'T))^ which corre- 
sponds to the microscopic density) and initial p.d.f. for the position of the particles /, the 
correlation between At^ and Bj-^ is defined as 

{Ar^Br.f = {Ar^Br,) - {ArJ{Br,) (2.3a) 

where 

i^Ta^n) = j^dx\\ . . j^dx'J f{xW . . . ^x'J) 

, . . . , Xj^ \xi , ■ ■ ■ , Xj^ ; ra) 

X ^ dx? ... 1^ B{{xf})G{x['\ . . . , 4) \X^^\ . . . , - Ta) 

(2.3b) 



{ArJ = j^dx^? ... j^dx%^ f{xf\...,x%^) j^dx^... j^dxNA{{x,}) 

G{xf', . . .,x^^\xi, . . .,XN;ra), (2.3c) 

and similarly the definition of {Br^). The quantity {Ar^B^-^) is referred to as the distribu- 
tion function. 

In the Langevin equation formalism the correlation between At at time r = and B^- 
at time is defined as 



{A^Bt,)^ = {AoBt,), - {Ao),{BtX (2.4) 

where ( )o denotes the average over the initial particle distribution, while the overline de- 
notes the average with respect to the random force for times from to r. The equivalence 
between the Langevin and Fokker- Planck formalisms implies 



{A^Bry = {AoBt,}, (2.5) 

When using the Fokker-Planck formalism we will make use of the general fact [11] 
that 



~ h^PdX] A^dXj) ^2 



pax] 4V9Aj.y '2 dX^ 

1 N 

-Enjn,. (2.6a) 

p j=l 



where 



Note that 



^ 1 d ipdW , 



n .g-w = 0. (2.7) 
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The formula (2.6) shows that after conjugation with e~^^^'^ the Fokker-Planck operator 
transforms into a Hermitian operator. For some particular W , which include (1.2), we 
have 

^pw/2Ce-^w/2 ^ _hjj _ (2.8) 

P 

where H is the Schrodinger operator for a quantum mechanical system with one and two 
body interactions only, and Eq is the corresponding ground state energy. For example, 
with W given by (1.2a) the equation (2.8) holds with 

^^ = -E£. + tS-?+W/2-1) E T^^-^r (2-9) 

i=i ^-^j ^ j=i i<j<k<N y-^j -^1^) 

Note from (2.6a) and (2.8) that the ground state is (up to normalization) given by 

^0 = (2.10) 



2.1 The current- current correlation 

In the Langevin description of Brownian dynamics the current-current correlation is de- 
fined by choosing the observables At and Bt in (2.4) as the classical current: 

A, = jo(^a), Bt = jV(x,), Ux) X: '^^5{x - x^{t)). (2.11) 

In the Fokker-Planck description the classical current has no immediate meaning as the 
velocities do not explicitly occur in (1.1). To define the current in this situation we make 
use of the formula (2.6a) and rewrite (2.3) in terms of time-dependent operators. 

For this purpose we note from (2.2) that 

JV 

G{xf\ ...,xP . . . , r) = e^^/^ U 6{x['^ - x^) (2.12) 

1=1 

where the Fokker-Planck operator acts on {xp"*} only. Substituting (2.12) in (2.3b) allows 
the integration over {.x['^''} to be carried out. Then substituting (2.12), with {a^l^-*}, {x;"^^} 
replaced by {a;^^''}, {a;^'*}, in the result and integrating over {a;^''} we obtain 

{AM = I dxf^ ...^ dx%^ 5.,e^(-^--»)/^A,„e^-»/V(a:?\ ■ ■ • , xg^). (2.13) 
Now it follows from (2.6a) that 

^Crh ^ g-/3M//2g-r^^''^^ntn^//37g/3H//2^ (2.14) 

so we can rewrite (2.13) as 

{AM = dxf^ ..-l dxP e-^^/2S(r,) A(r„)/e^^/^ (2.15a) 
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where 

A{r) := e^E,^intn,//9.^^^-.E,1.ntn,//9. (2.156) 

and similarly the definition of B{t) (in deriving (2.15a) we have used the fact that 
e-Ef=in]n,//37g-/3w^/2 _ ^^-^^ fou^^g ^^^^ ^2.7)). 

Using (2.15b) we can define the current jrix) by the continuity equation 

^n{x,T) ^ —^j{x,r), (2.16) 
where n{x,T) is defined by the r.h.s. of (2.15b) with 

N 

Ar = ririx) := J2 ^(^ - (2-17) 

i=i 

which is the microscopic density. Now (2.15b) is the imaginary time quantum mechanical 
formula for the evolution of the operator Ar in the Heisenberg picture. Further the 
Hamiltonian J2jLi njllj is of the form 

-EI^ + '^(^1'---'^^)' (2-18) 

3=1 

SO we conclude from the usual quantum mechanical verification of the continuity equation 
that the sought definition of the current is 

M-) - -^E {-^g^/i- - -Ar)) + Six - xAr)).g^). (2.19) 

As written (2.19) assumes the particles are on a line. However the same formula applies 
for each component of the current vector in higher-dimensions if we simply replace jrix) 
by j^"^(x) and d/dxj by d/dx^"\ 

Note that with / proportional to e~^^ (the equilibrium state) in (2.15a) we can write 

{AM = {HB{n)A{Ta)\iJo) (2.20) 

which is precisely the quantum mechanical formula (with r = itf3'y) for ground state cor- 
relations. This shows that, up to a factor of (i/7/3)^, the current- current distribution for 
Brownian motion described by the Fokker-Planck equation (2.1) with initial p.d.f. given 
by the equilibrium distribution is identical to the current-current distribution of the corre- 
sponding quantum mechanical system (recall (2.8)). The equivalence between the density- 
density distributions in this situation has previously been shown by Beenakker and Rejaei 

[!]• 

2.2 /-sum rule 

Denote the density-density correlation by S so that 

N N 

SiiXa, 0), {Xb, n)) = ( E SiXa - XjiO)) E ^i^b - Xk{n))) (2.21) 

j=l k=l 
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where the average on the r.h.s. is defined according to (2.3). When the initial distribution 
/ in (2.2) is proportional to the equilibrium distribution e"^^ we have just remarked 
that S is identical (with r = it(3^) to the ground state density-density correlation of 
the quantum system with Hamiltonian J2jLi n^H^. Since this Hamiltonian is of the form 
(2.18) it follows that S must satisfy the so called /-sum rule [10]: 

= -^P (2.22a) 

:0 7P 

where ^ 

S{k, r) := / Siixa, 0), {x^, r))e"'(-'^--»^ d{xa - x^) (2.22&) 

and it is assumed that the ground state is translationally invariant so that S depends 
on Xa — xi,. Subject to this assumption (2.22a) is valid for N finite as well as in the 
thermodynamic limit. 

In fact it is possible to derive and indeed generalize the /-sum rule entirely within the 
Brownian motion setting using either the Fokkcr-Planck or Langevin equation formalism 
(for definiteness we will consider the latter). The generalization is that we can take the 
initial distribution as proportional to the Boltzmann factor e'^ ^ (where we may have 
(5' ^ P) and we do not assume translational invariance. 

We consider the quantity 

S{x, k, r) := / Siix, 0), (x,, r))e''^» dx,, (2.23) 

— oo 

where S is given by (2.20) with the average now computed according to the r.h.s. of (2.5) 
because we will use the Langevin formalism. Differentiating with respect to r gives 



1=1 



(2.24) 



As T approaches zero, we see from the Langevin equation (2.1) that 

dW 



-fXj(r) = 



dxi 



T=0 



+ ^,(r) + 0(r) 



and so 



^ikxi{r) _ gifcx((0) / 



(2.25a) 



(2.256) 



where the term proportional to r not explicitly written is independent of ^jir). Sub- 
stituting (2.25) in (2.24) and using (2.1b) to compute the average with respect to the 
random force gives 



7- 



dS{x, k, t) 



dr 



= E^(5(x-x,(0))e^'=-(-)(-|^L^^+^ l'sir-r')dr))^ (2.26) 



The delta function in (2.26) is non-zero on the boundary of the interval of integration so 
we take one half of the value of the integrand. After averaging over the initial distribution, 
the corresponding term of (2.26) is then identified as proportional to S{x, k, 0). To simplify 
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the term involving dW/dxj in (2.26) we recall that the initial distribution is assumed 
proportional to e"^^ . Since 



dxj P' dxj 

integration by parts gives that this term can be rewritten as 

-f^(^>^>0) + l^(e^'=V(x)). (2.27) 

Combining these results gives for our generalization of the /-sum rule 

^^l^ = -^1^ - 0) . l^^(e-pW). (2.28) 



which reduces to (2.22a) in the case P' — (3, p{x) — p. 

We emphasize that (2.28) holds for a finite system, as well as in the thermodynamic 
limit. It has been presented in the one-dimensional case, but by introducing the higher- 
dimensional Fourier transform in (2.23) and repeating the working we see that (2.28) 
holds with k"^ replaced by Y,a{k^°'^T the first term and k, d/dx, kx replaced by k'^°'\ 
d/dx^'^\ k^'^^x^"^ (with summation over a) in the last term. 



2.3 Static current- current distribution 

For the particular Schrodinger operator 

«--|:g^W/2-l)(f)\j^^-.^^J^. 0<.,<i. (2.29) 

which is related via equation (2.8) to the Fokker-Planck operator (1.1) with potential 
(1.2b), it has been shown by Taniguchi et al. [3] that the static current-current distribution 
is given in terms of the static density-density distribution by 

l(lpQ\j{Xb,Q)j{Xa,^)\ijQ) = ■ 2 ^^^^ . /r {'^o\p{Xb,0)p{Xa,0)\'llJo), Xa 7^ Xb- (2.30) 

sm Tr{xb - Xa)/L 

From the remarks below (2.20), (2.30) is equivalent to the statement that for the Fokker- 
Planck system with potential (1.2b) and initial p.d.f. equal to the equilibrium distribution, 

. . , , . , jdin/LY / ^ \2/ / X / x> /„ 

{jo[Xa)jQ[Xb)) ^ T^— VTt[—5) {'>T'o{Xb)no[Xa)) . (2.31) 



sin TT{xb - Xa)/L ^7/? 

The equation (2.30) was derived in ref. [3] using the factorization (2.6a). The same method 
readily extends to provide an analogous result for all Fokker-Planck systems in which W 
consists of one and two body potentials only. 

Thus we are considering the average (2.15a) with Ta and T(, equal to zero, the observ- 
ables A and B given by the current (2.19) and the initial distribution / proportional to 
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e To apply the method of ref. [3], we substitute for d/dxj in (2.19) using Ilj as 
defined by (2.6b) to obtain 



io(^) = E {^W - ^j) + - (2-32) 

PI j=i 



Prom this formula and (2.7) we have 
where 



{jo{xa)jo{x,)) = (-f )' E i / d^?^ ■■■I dx^N e-^^/^<5(x, - Xj)njUlS{xa - x,)e-^^/2 

(2.33) 

Let us now suppose Xa 7^ X},. Then the j — k term vanishes. But for j ^ k we have 

Using this formula, the fact that 11^5(3; — Xk) = S{x — Xfe)!!^ for j ^ k, the property (2.7) 
and the definition (2.17) we deduce from (2.33) that 



^1— 



{no{xa)no{xb)), (2.35) 



which is the sought formula. We remark that it is also possible to derive (2.35) within the 
Langevin equation setting used in the derivation of the /-sum rule given in the previous 
section. Also, for W given by (1.2b), note that (2.35) reduces to (2.31). 

Our derivation has been presented in the one-dimensional case. The same result 
applies in higher dimensions for the correlations between components jjf\xa), jif\xi)) of 
the current vector provided the Xi,Xa and X2,Xb outside the average on the r.h.s. are the 
components (a) and (/5) respectively of the corresponding position vector. 

2.4 Small- A; behaviour of structure factor 

In this section we will consider Brownian motion described by the Fokker-Planck equation 
(1.1), with the potential W consisting of general one and two body terms 

N 

W = Y.V^{xj)+ Y: V2{\xj-Xk\) (2.36) 

i=l l<j<k<N 

such that the potential V2 is non-integrable at infinity. For the particular case of this type 

V^ix) ^ -log\x\ (2.37) 
Dyson [5] introduced the hydrodynamic equation 

7j(x,r) = -p{x,t) — {v^{x) - J_^dx'p{x',r)V2{\x-x'\)) (2.38) 
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(note that the r.h.s. represents the force density) to study the large-wavelength density 
fluctuations of the system. Here j{x,T) and p{x,t) represent smoothed, continuum ap- 
proximations to the microscopic current and density. 

In the case (2.37) it has been deduced in ref. [1] that for a uniform initial state 

p{x',0)=p (2.39) 

the Fourier transform of the density-density correlation deflned by (2.22) has the small-A; 
behaviour 

S{k,T) ~ 5(A;, 0)e-*^'^''^^('=)/^ (2.40) 

where ^ 

V^(k) := r ¥2(1x1)6'"^ dx. (2.41) 



Inspection of the derivation given in [1] shows (2.40) to be a consequence of (2.38) inde- 
pendent of the particular V2, so its validity for general V2 is tied to the asymptotic validity 
of the hydrodynamic approximation (2.38). We expect (2.38) to correctly describe the 
large- wavelength density fluctuations whenever the potential V2(|a;|) is not integrable at 
inflnity (for integrable potentials the 'electric fleld term' involving the integral on the 
r.h.s. of (2.38) is not present, rather the force density is due to a pressure gradient; see 
the next subsection). Now for potentials not integrable at infinity it is expected [13] that 

S(k,0) i (2.42) 

^ ' mk) ^ ' 

(here (5' is used as in the discussion of the generalized /-sum rule, and thus may be 
different to (3 which occurs in (1.1)) and so (2.40) reads 

Sik.r) i g-fcW2(fc)/7_ (2.43) 

Note that this expression is consistent with the generalized /-sum rule (2.28) with p{x) = p 
(note that the first term in (2.28) is of a lower order than the second since, as is seen from 
(2.42), S{k, 0) — > as /c — >• whenever V2(|x|) is not integrable at infinity). It is expected 
to be valid for non-integrable potentials in higher dimensions provided (2.41) is suitably 
modified. 

One consequence of (2.43) is the sum rule 

/ S(k,T)dT 1 as A;^0. (2.44) 

Jo ^ ' ^ pl3'k\V2{k)Y ^ ' 

In the particular case when V2 is given by (2.37) and /3' = /3, we noted in Section 2.1 
that S(k^ t) in the Fokkcr-Planck system is identical to S{k, r) for the ground state of 
the Calogero-Sutherland quantum system with Schrodinger operator (2.9). Now for a 
quantum fiuid the integral (2.44) is related to the compressiblity of the ground state. 
This is consistent with the fact that for V2 given by (2.37) 



Mk) ~ ^ (2.45) 



and so the r.h.s. of (2.44) is a constant. 
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Another consequence of (2.43) is the evaluation of the large-r mean square displace- 
ment of a particle from its initial position. Thus from ref . [14] , assuming this displacement 
diverges, S{k,0) diverges at the origin and that the state is homogeneous, we have for 
large-T 

{{x{t) - - ^ X i ("^(^^ 0) - Si^^ ^)) (2-46) 

where c <^ 1 is some constant. To make further progress, suppose V2(r) ~ r~", < a < 1 
as r ^ oo, so that [15] 



.la-l 



V2{k) - -r as ^ 0. (2.47) 

T[a) cos 7rQ;/2 

The bounds on a are required for the validity of (2.46). Indeed the bound a < 1 is 
required so that S{k, 0) diverges at the origin, while substituting (2.47) in (2.43) and then 
substituting the result in (2.46) and shows that for (2.46) to diverge as r — > oo we must 
have q; > 0. To calculate the large-r behaviour we make the substitutions of the previous 
sentence and change variables y — k^^^ii: /V{a) cos7ra/2^rp/^. This gives 

((a;(r) -a;(0))2) ~ c(/5',7,p)r"/("+') as r ^ oo (2.48a) 

where 

. ^(PZ2)!;p( r(«)cosW2 )V..... r _ 

P'ttp^ a + 1 ^ TT ^ Jo 

(2.486) 

which is the sought formula. In the case a > 1 (integrable potential) it is shown in ref. [14] 
that {{x{t) — x{0)y) is proportional to r-*^/^, independent of a. 

For a one-dimensional translationally invariant initial state which is the equilibrium 
state, the small- /c behaviour (2.43) also implies the small-/c behaviour of the Fourier 
transformed current-current correlation 



C(^,t):= / {jo{0)jr{x)fe"^^dx. 

J —oo 



To sec this, we recall [1] that for a one-dimensional translationally invariant system, a 
consequence of the continuity equation (2.16) is that 



where 



/'■OC 
S{{0,Ta),{x,n))e^'^^dx. 
-oo 



But by the assumption that the initial state is the equilibrium state S{k; Ta, Tf,) — S{k, Tb — 
To) and so (2.49) reads 

Substituting (2.50) in (2.43) (with = (3) gives for the small-fc asymptotic formula 

C{k,T) -^A;2^(fc)e-*^'^''^2«/^. (2.51) 

10 



For the Dyson log-gas, when V2{k) is given by (2.45), (2.51) gives the r = result 

'5(M, ^ 

By taking the inverse transform we deduce the large-Xab behaviour 

which is consistent with the large-Xa^ behaviour of the thermodynamic form of (2.31) , 
since in this limit (?^o(3^6)^o(^a)) ~ ■ 

We also have the analogue of (2.44) in the setting of the validity of (2.50). Thus 
integrating (2.50) with respect to r gives 

where to obtain the last equahty the /-sum rule (2.22a) has been used. Note that unlike 
(2.44), the equation (2.53) is valid for all values of k. 



2.5 Second order correction 

To obtain the next order correction in the small-/c expansion (2.40) we refine the hydrody- 
namic equation (2.38) to include a pressure gradient —dp{x,T)/dx as an additional force 
density on the r.h.s.. The strategy for the solution of this equation is the same as given 
in [1] for the solution of (2.40). The first step is to differentiate both sides of the equa- 
tion with respect x and to then substitute for dj{x,T)/dx using the continuity equation 
(2.16). After linearizing in Sp{x,T) := p{x,t) — p (note that the pressure is linearized 
by p{x,t) — p = {dp/dp)Sp{x,T) where p denotes the bulk pressure) and assuming the 
equilibrium condition that for large-L 

d 



— (v,{x)-pJ^V2{\x-x'\)dx')^0 



we obtain 

dSp{x,T) p d"^ /"°° T^n / / \ J / ^dp d"^ 



J V2{\x - x'\)6p{x', t) dx + ~^^^^(^' ^) 



dr 7 dx'^ J-oo "f dp dx"^ 

Taking the Fourier transform of both sides and integrating by parts twice gives a simple 
differential equation for 5p{k, r) which has solution 

8p{k, t) = 6p{k, 0) exp ( - (k^pV2{k) + ^) . (2.54) 

But 

Sik,T) = j{6pi-k,0)Spik,T)). (2.55) 
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Substituting (2.54) gives the desired second order correction to (2.40): 

S{k,r) ^ Sik,0)exp(^~ (k" pV2 {k) + ^k"") ^) . (2.56) 

Note that for consistency with the second term of the /-sum rule (2.28) with p{x) — p, 
which gives the leading small- /c behaviour, we must have 

T-~ — — VTTX- (2.57) 

A thermodynamic derivation of this result is given in the appendix. 



3 Exact results 



3.1 GOE GUE transition 



Consider the Dyson Brownian motion model specified by the Fokker-Planck equation 
(1.1) with P — 2, W given by (1.2a) and the initial distribution proportional to . 
Since is proportional to the eigenvalue p.d.f. for the Gaussian Orthogonal Ensemble 
(GOE) of random real symmetric matrices, while e^^''^ is proportional to the eigenvalue 
p.d.f. for the Gaussian Unitary Ensemble (GUE) of random Hcrmitian matrices, this 
specific Dyson Brownian motion model describes the GOE GUE transition. In the 
N ^ oo limit (suitably scaled) the system is translationally invariant with bulk density 
p and the Fourier transform of the density-density correlation (2.22) is given exactly by 
[4, eq.(3.32) with the correction that an additional factor e(^^l'=lp/T)(i-l'=l/2'^^) be included 
in the first term of (3.32c)] 



< |A;| < 2np (3.1a) 



27 

S{k,T) = — —T exp(— r7r| A;|p/7)sinh(rA; 727) 
7rT|A;| 



l-l+\k\/np I 



27r Jl ki 



27 

S{k,T) — — -— exp(— TA;^/27)sinh(T7r|fc|p/7) 

7lT\k\ 

M^rkV2, r^+l'l/^^ dh^^-^^mk.h^ \^ > 27rp (3.1b) 
27r J_i+|fc|/^p ki 

This result provides an illustration of the generalized /-sum rule (2.20) in the case (3 7^ j3' . 
Thus according to (2.28) with p{x) = p, (3 = 2, (3' = 1 and x = 0, for all k (3.1) must 
satisfy 



dS{k,T) 



an ^ ^^P 



OT r=0 ^7 7 

An elementary calculation verifies that this is indeed so. Also, the small-A; expansion of 
(3.1a) is consistent with (2.56) and (2.57) with (3' = 1, V2{k) given by (2.45) and dp/dp 
calculated from the equation of state (3'p = (1 — j3'/2)p. 
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3.2 Free fermion type system 



As remarked earlier, for a number of special choices of W including (1.2), the equation 
(2.8) holds relating the Fokker-Planck operator to a Schrodinger operator. A further 
property of this Schrodinger operator, as is illustrated by (2.9), is that at the special 
coupling /3 = 2 the coefficient of the two-body term vanishes and the system can be 
regarded as free fermions in an external potential. 

Consider such a situation in which the one-particle Schrodinger operator. Hi say, 
has a complete set of real, orthonormal eigenfunctions {ipk}k=Q,i,... and corresponding 
eigenvalues {e^ -|- £'o}fe=o,i,.-- the particles are confined to an interval /. A standard 
calculation (see e.g. ref.[l]) gives that the density-density correlation for the A'"-particle 
system is given by 

oo iV-1 

5((x,o),(x',r)) = J2M^)M^>~''''^'''T.M^)M^'y'^^''' 

p=0 q=0 

N-l N-l 

p=0 q=0 

Since here we have set t = r/i^j, this represents the density-density correlation for the 
corresponding Fokker-Planck system with initial state equal to the equilibrium state. This 
exact result can be used to illustrate the generalized /-sum rule (2.20) in the case of finite 
N and non-uniform density p{x). 

Thus we want to show that (3.2) satisfies 



dr 



where S{x,k,T) is defined by (2.23) with S therein given by (3.2), and p{x) is given by 

N-l 2 

P(^) = E (M^)) (3-4) 

p=0 

(this latter formula follows from the free fermion calculation). For this purpose we note 
that differentiating the second term in (3.2) with respect to r and then setting r — gives 
zero as the two terms which result cancel. Using this feature and the fact {Hi—EQ)ipq{x') — 
eqipg{x') we see that 

a 1 oc N-l 

-S{{x,0),{x',t))1^^ = -— EV'p(^)((^^i- Wp(^')) E^.W.(^') 

1 oo N-l 

+ -5J2M^)M^') E M^)iHi - Eo)i^g{x'). (3.5) 

'^y p=0 q=0 

Taking the Fourier transform of both sides with respect to x' and integrating by parts 
(twice) the first term on the r.h.s. using the explicit formula 
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we see that the last term cancels. Furthermore, noting by completeness that 



J2^p{x')'tpp{x) = 6{x' - x) 

p=0 



we obtain 



^S{{x,0),{x',r)) 



T=0 



q=0 
N-1 



-—^ I - x') Y: dx' 



(3.6) 



where to obtain the last equality (3.4) has been used. Comparing (3.6) with the /-sum rule 
formula (3.3) shows that the required agreement with that result has been demonstrated. 



3.3 The Dyson Brownian motion model for rational (3 

In the thermodynamic limit the density- density correlation S{x. t) for the Dyson Brown- 
ian motion model specified by the Fokker-Planck equation (1.1) with W given by (1.2b) 
has been calculated for all rational values of (3. With (3 /2 :— \ — p/ q [p and q relatively 
prime) the result is [9] 

lim Six, t) 

JV,Z,-.oo ^ ' 
N/L=p 

= CpM n / dxi^ dyjQl^qF{q,p, X\{xi,yj}) cosQp^qX exp{-Ep^qT/2X-f) 

i=l "^0 j=l "^0 

(3.7a) 

where the momentum Q and the energy E variables are given by 

q p q p 

Qp,, - 27rp(^a;, + ^p,</ ■= C^^pfi^Z^pi^j) + J2^H{yj)) (3.76) 

1=1 j=l i=l j=l 

with 

ep{x)^x{x + X) and eniv) ^ Xy{l - y), (3.7c) 
the form factor F is given by 

and the normalization is given by 

^ _ A^^fa-^)r^(p) r^(A)rni/A) 

"^^■^^-"^ 2n^p\q\ YIU T\p - X{i - 1)) m=i T\l - U - 1)/A) • ^^-'"^ 
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This exact result can be used to illustrate the sum rules (2.43) and (2.57). To begin, 
note that S{x, r) is even in x and so S{k, r) is even in k. It therefore suffices to consider 
the case A; > 0. With this assumption, from (3.7) 

<i POD ^ /■! 

= nCp,g{X) Yl dxiYl dyjQl ,jF{q,p, X\{xi, yj}) 5{k - Qp^q) ex.p{-Ep^qT/2X-f). 

Changing variables Xi ^ kxi, yj i— > kyj gives 
S{k,r) ^ 7rkCp,g{X)e-''''^/^ 

X ]J / dxij{ dyjQl^gF{q,p, X\{xi, yj}; k) 5{l - Qp,q) exp(-Ep,g,fcT/2A7). 

(3.8a) 

where 

F(,, p, A| {X., y,}; ^) = R ft (x. + Xy.^ ^^<^'\^^- -£^^<r\y- Vr?'' 

i=ij=i nLi [xi{kxi + A)) UU {Xyj{l - yjk)) 

(3.86) 

and 

Ep,g,, = (27rp)V(^^- - AE^D- (3-8c) 

i=i j=i 

In the limit A; — > the integral in (3.8a) is independent of k, so to leading order we 
have 

S{k,T) ~ 7r|A;|Cp,,(A)/(A)e-l'^l'^'-/^ (3.9) 

where 

POO P POO 

/(A) = n/ dxiH dyjQlqG{q,p,X\{x„yj})6{l-Qp,g) (3.10) 

i=l "'O j=l ■'^ 



with 

q p 



G(g,p,A|Ky,}) = U Ei-^ + Ay.)-^ "-"' "Z"^; '^va'^''''' - (^-H) 



To evaluate (3.10) we rewrite it as 

1 



;>oo y /-oo POO 

/(A)= lim — / (iwe'" [] / 11 / c?%<3L<^(g,P, A|{a;i, |/j})( 

(3.12) 

Now change variables Xi i— > j^Xi, yj i— > j^yj- This shows 



(3.13) 

The first integral above equals unity in the limit e — 0+, while the multiple integral has 
been evaluated in an earlier work [16] as equal to l/(Cp,q(A)7r^/3). Substituting in (3.13) 
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and then substituting the result in (3.9) gives the hydrodynamic result (2.43) with the 
substitution (2.45). 

It is also possible to illustrate (2.57) by expanding (3.8a) to the next order in k: 
F(q,p,X\{xi,yj};k) ~ G(q,p,X\{xi,yj})(l + k{l - l/X)Qpj27rpy 

This gives 



S{k,0) 



\k[ 



+— (l-l/A)7rA;2Cp,,(A)ny^ IT ^ dyjQl^G{q,p, X\{xi,yj})6{l - Q^,,). 

(3.14) 



Comparing the multiple integral above with that in (3.10) shows that they arc identical 
apart the factor of Q^^ which reads Qpg in (3.10). But this factor docs not effect the 
value of the integral due to factor of 5(1 — Qp,q). Thus from the above working we have 

This is precisely the result (2.57) expanded to second order in k with /3' = /3, the substi- 
tution (2.45) and the partial derivative evaluated from the equilibrium equation of state 



4 Conclusion 

We have presented a systematic study of exact and asymptotic formulas relating to dy- 
namical correlations for ovcrdamped Brownian motion. Our main new results are the 
operator formula (2.19) for the current in the Fokker-Planck description, the generalized 
/-sum rule (2.28), the formula (2.35) relating the static current-current correlation to 
the static density-density correlation and the asymptotic formula (2.52) for the small-A; 
behaviour of the dynamical structure factor. We have illustrated these formulas using 
known exact evaluations of the dynamical density-density in certain special models. 

Our work builds on the results presented in refs. [1-3]. In this regard we point out 
that in ref. [1] the current is defined not in terms of dxj{T)/dT as in (2.11), but rather 
with this factor replaced by dxj{T) / dr^^"^ . The latter choice is appropriate in applications 
of the Dyson Brownian motion model to chaotic spectra as the perturbing parameter X 
is related to r by r = X^. 
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Appendix 



Here we will derive (2.57) using a thermodynamic argument. For notational convenience 
we will assume the domain is an interval of length L, although our argument applies 
equally well in higher dimensions. The setting is a uniform one-component system with a 
pair potential non-integrable at infinity. Such a system requires a neutralizing background 
for thermodynamic stability, and supposing the particles carry unit charge, the particle 
and charge density fluctuations are identical. We consider the microscopic density and 
suppose it has a fluctuation of the form 



5p{x) = pfce''^^ + pie 



—ikx 



(Al) 



Here pk is 1/L times the Fourier transform of the microscopic density. If 2T[/\k\ is a 
macroscopic wavelength, p^ is a good collective variable which can be used instead of the 
particle coordinates. The system will gain an electrostatic energy 



1 (■L/2 i-L/2 

- / dx'5pix') / dx5pix)V2i\x-x'\). 

2 J-L/2 J-L/2 



(A2) 



Let f{p) be the free energy per unit length of the locally neutral system. To second 
order in 8p{x) the density fluctuation changes this free energy by an amount 



dp p=po 2 op'' 



P=PO 



(5p{x)y. 



(A3) 



Hence the total change in free energy due to the fluctuation is 

1 1-^2 , , _ /-VS , 1 



L/2 



p=Po J-L/2 



dx (6p{x 



P=PO^ 



(A4) 



In general the probability density function for an event in the system is proportional 
to where 5F is the corresponding change in the free energy. From (2.55), where 

now 6p{k) — Lpk, we have 

m = L{\p,\'). (A5) 
Computing the r.h.s. of (A5) from the Gaussian distribution e~^'^^ gives (2.57). 
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